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ABSTRACT 

Based on a distributed parameter model for vibrations, an approximate finite 
dimensional dynamic compensator is designed to suppress vibrations (multiple modes 
with a broad band of frequencies) of a circular plate with Kelvin-Voigt damping and 
clamped boundary conditions. The control is realized via piezoceramic patches bonded 
to the plate and is calculated from information available from several pointwise observed 
state variables. Examples from computational studies as well as use in laboratory 
experiments are presented to demonstrate the effectiveness of this design. 
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1. INTRODUCTION 


In recent years, a great deal of research has been carried out on the development 
and derivation of control designs from an infinite dimensional state-space approach; 
however we are unaware of any implementation based on such designs being reported in 
the literature. We have presented briefly in [1] some preliminary experimental results on 
implementation of an output feedback control which was designed based on an infinite 
dimensional (or distributed) system. In this paper we attempt to summarize the design 
methodology and give further discussion of the implementation. 

The feedback control system was implemented on a circular plate with a piezo- 
ceramic patch as actuator. This choice of structure was motivated by the fact that it 
is an isolated component from the structural acoustic system described in [7, 8]. I he 
structure in that system is made up of a hard walled cylinder with a (damped circular 
plate at one end and the control problem consists of using piezoceramic patches on 
the plate to reduce t he interior structure-born sound pressure levels which result when 
the plate is subjected to a strong exterior acoustic field. The partial differential equa- 
tion (PDE) system which describes the dynamics of this circular plate is presented in 
Section 8 below. 

In our control design, three primary concerns are: 1) presence of disturbance 1 in 
both input, and output of the system; 2) robustness of control: 3) lack of full state 
measurement. Those concerns lead to a design problem involving dynamic compen- 
sators for distributed parameter systems. A great deal of recent research has been 
carried on the individual or combined topics of our concerns here, for example, see 
[2. •] , 13, 15, 1(), 17, 18, 20, 21], and the references therein. This problem involves difficult 
issues and many theoretical and computational questions remain to be resolved. 1 he 
purpose of this paper is to demonstrate how finite dimensional control theory together 
with approximation theory for certain optimal control problems can be used to success- 
fully design and implement feedback controllers for flexible struct ures. A finite dimen- 
sional dynamic compensator design is outlined in Section 2. The approximation scheme 
which leads to a finite dimensional control problem will be discussed in Section 1. A nu- 
merical example is given in Section 5 along with experimental results to provide 4 prelim- 
inary validation regarding the implementation of a PDE- based (distributed parameter 
or infinite dimensional system based) method for reducing structural vibrations. 

2. CONTROL PROBLEM FORMULATION 

We first consider an n- dimensional system 

2/(0 = Ay{t) + Bu{t), y{ 0) = y 0 , 

!U(t) = Cy(t ), (1) 

z(t) = Hy(t) + C;u(t). 

where the state variable y is in //T\ the control u is in IR V \ the measurement y t)h is in W 
and the controlled output c is in Iff for some finite positive integers n, m, ]h and ?*. 
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The coefficients A , B, C, H , and G are time invariant matrices. The performance 
index (or cost function) is given by 


J(u) 




(It 


J Q |(W0iy(0) + (Ru{t),u(t))^dt, 


( 2 ) 


subject to (1). In (2), Q = //'//, R = G*G , where we assume H f G = 0. The control 
problem is to find a controller u £ oo; iR m ) which minimizes the cost function (2). 

We are interested in the case when, as in most practical situations, measurement 
of the full state is not available (p < n). (We note that when y 0 b{t) = y{t) the solution 
can be obtained by applying the well known linear quadratic regulator (LQR) optimal 
control theory.) One possible approach is to build a state estimator or observer to 
reconstruct the state from the measured partial state. One can then feed back this 
reconstructed state. In this paper, we consider a full order observer for simplicity (it is 
adequate for the plate experiments described below). Design of low order observers is 
important when the dimension of the control system is large, and we refer to [10, 15] 
for details on reduced order observers. 

Let the reconstructed state be denoted by y c (t). We consider the standard com- 
pensator (of Luenberger type [19]) for system (1) given by 


y c (t.) = A e y c (t) + F y ob (t), 

A c = A- FC - BK , (3) 

u{t) = -Ky c (t), 

for a properly chosen feedback gain K and observer gain F so that the reconstruction 
error | y c {t) — y{t ) | — » 0 as t — > oo and the closed loop system 


y{t) 


A -BK ' 

y(t) 

Vc(t) _ 


FC A c 

Vc(t) 


is exponentially stable. Intuitively, we would like to choose F such that the observer 
poles of A — FC are deep in the left half complex plane to obtain fast convergence of 
the reconstruction error. This must be done with care since an observer so constructed 
is very sensitive to any observation noise that may exist. 

Among several compensator designs, we first consider the so called optimal com- 
pensator. Suppose that the matrix Q is nonnegative-definite, R. is positive-definite, the 
pair (A y B) is stabilizable, (A, (7) is detectable, (A,G) is controllable, and (A,//) is 
observable. Then there exist unique (minimal) optimal feedback gain K and observer 
gain F given by 


= r-'b'u 

- PC'R- 1 


K 

F 


2 


( 4 ) 

( 5 ) 



where n and P are unique nonnegative-definite solutions to the following regulator and 
observer algebraic Riccati equations 

n A + A'U - n BR-' B'n + Q = 0, (6) 

PA 9 + AP- PC PC' CP + Q = 0, (7) 

respectively. Thus the optimal estimator is obtained and given by (3)-(5). In (7), Q is 
a nonnegative symmetric matrix and R is a positive symmetric matrix. The matrices 
Q , /?, Q, and R are determined by some design criteria for the specific control problems. 
We point out that this “optimal” observer can be defined without depending on the 
(traditional) stochastic formulation. The name “optimal” is derived from the stochastic 
interpretation of the above design (see [13, 19] for further discussions). Briefly, the 
above described observer 

2 Jc(t) = (a~ BRT'B'Ujydt) + PC’fr l c(y(t) - y c {t )) (8) 

is simply the Kalman- Buc.y filter if we consider the system (1) disturbed by the uncor- 
related stationary Gaussian white noise V\ (t) and v 2 (t): 

y(t) = Ay(t) + Bu(t) + v-i(t), ?/(0) = y 0 , 

y 0 b(t ) = Cy(t) + v 2 (t), (9) 

z(t) = Hy(t) + Gu(t), 


where 


= 0, £{v\{t) u((r)} = QS(t. - T ) 

£{v2(t)} = 0, £{v 2 {t)v' 2 (T)} = RS(t - t). 

Here £{ } is the expected value. The observer (8) is optimal in the sense that the limit 
of the mean square reconstructed error 

Km £{(y(0 - y e (tj) w (j /(<) - 2 /c( 0 )} 

( W is a weighting matrix) is minimal with respect to all other observers (e.g., see [19]). 

Even though the optimal compensator provides us with the desired performance, 
it is well-known that it may lack robustness. To design a robust dynamic compensator, 
let us consider the system (1) with input and output disturbance iv(t) 

y(t) = Ay(t) + Bu{t ) + Dw(t), y( 0) = y 0 , 
yM = Cy(t) + Ew(t ), (10) 

z(t) = Hy(t) + Gu(t.), 

where the disturbance vector rw(^) is in IR q for some finite positive integer q. The 
coefficients D and E are time invariant matrices. Furthermore, we will restrict ourselves 
to matrices //, G, D, and E such that H' H — Q > 0, G'G — R > 0, H'G - 0, DO' — 
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Q > 0, EE ' — R > 0, DE f — 0. The more general ease where the cross product 
terms H'G and DE f are not zero can be dealt with in a similar manner with slight 
modifications (see [9]). 

Our objective is to design a robust controller that provides acceptable performance 
with disturbed incomplete state measurements. One such design technique is the so- 
called Hoo/MinMax compensator given in [9]. One formulates -control problems 
in the time domain and obtains a soft-constrained dynamic game associated with the 
disturbance attenuation problem. The control problem is formulated as a form of 
optimization of a performance index (or cost function). For this purpose, we introduce 
the extended performance index: 

= jf (|jr(0| J ~ 7 2 \ W U )\ 2 } dt 

° , ( 11 ) 

= J o + (Ru(t),u(t)) - 7 2 ( u .’(/,), iw (/))| di 

subject to ( 10). The optimization problem is to find a controller u* £ U = L-2(0, oo;ZR m ) 
and disturbance iv* £ W = L*2(0,oo; IR (] ) such that 

J* = inf sup J 7 (u, ic) = J 7 (uVuF). 

u ^ lJ w£W 

One seeks necessary and sufficient conditions on 7 so that quantity J * is finite. The 
lower bound of 7 for which J* < 00 is the optimal minimax attenuation level and is 
denoted by 7*, i.e. 

7* = inf {7 : J * < 00}. 

The first part of the optimization problem formulates the soft-constrained game and 
the second is a disturbance attenuation problem (7 is the attenuation level). It can be 
shown that the results of this optimization problem yields a bound for the //00-norm 
of the transfer function from disturbance w(t) to the controlled output -(/). 

To be more precise, the central results for this control problem can be summarized 
as following. Let the pair ( A , B ) be stabilizable, (A, C) be detectable, (A, G) be control- 
lable, and (/!, //) be observable. For a given attenuation 7 > 0, there exist (minimal) 
positive definite solutions FI and P to the following two algebraic Riccati equations 


UA + A'n - U(BR- { B f - 7" 2 Q)n + Q = 0, (12) 

PA! + AP - P{CfC'C - 7 ~ 2 Q)P + Q = 0, (13) 

respectively. Moreover, if the spectral radius p of PY 1 satisfies the condition 

p(PYl) < 7 2 , or 11 -7 2 /'- 1 < 0, (14) 

then there exists a unique optimal controller 

u(t) = -R- l B'Hy e (t), (15) 

and the state estimator y c {t) G Ei n satisfies 
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( 16 ) 


y c {t) = A c y c (t) + Fy ob {t), 

2/r(0) = y CO i 


where 



A e = 

A ■ 


F = 

(I 

In addition, we have 7 > 

7*- 


The resulting closed-i 

loop system 


2/(0 

' 


Veit) , 



K - FC + 7“ 2 Q1 

2 pm-' pr" o - 1 


A 

FC 


-BK 

A c 


?/(0 

Vc{t) 


with the controlled output 


:(t) = 


' H O' 


2/(0 

0 —GK 


. 2 / c (0 . 


(17) 


(18) 


is stable. 

Finally, if we let £(,s) and w(s) denote the Laplace transform of z(t) and w(t) 
respectively, then the transfer function from the disturbance w(t) to the controlled 
output z(t) is expressed by 


T( S ) = 


w(s) 


H 

0 


0 

-GK 


si - 


A -BK 
FC A c 




D 

FE 


(19) 


Furthermore, the norm of the transfer function (19) is bounded by 

nn-)iioo<7- 

Thus, if we follow this procedure we obtain a dynamic compensator which not only 
stabilizes the system with imperfect state measurements, but also provides robustness. 


3. STRUCTURAL MODEL 

In this section, the mathematical model used to describe the experimental setup 
is given. To reduce computational complexity, in our initial experiments the structure 
is axisymmetrically configured. We point out that all of the results and techniques 
presented here can be extended directly to the more general case of nonaxisymmet- 
ric configurations. The structure under study is a fixed-edge circular plate with a 
centrally placed circular shaped piezoelectric ceramic patch for actuation and sens- 
ing. The equations of motion will be formulated in polar coordinates (r, 0). Under 
the Love- Kirchh off plate theory with Kelvin-Voigt (or strain rate) damping, the trans- 
verse vibrations w(t , r, 0 ) of a plate of radius a subject to an axisymmetric external force 



g(t,r,6) are described by the system 


. n ,d 2 w dw d 2 M r 2d M r 

p{r ’ e) ~dF +c ' , 'm + lh :r + rli7 
w(t , a, 6) = 0, a, 0) = 0, 


1 dMe 

r dr 


= V 2 M pe + g 


0 < r < a 
t > 0 


( 20 ) 


where the internal bending moments are 


M r 

= D(r,d) | 

( d 2 w i/(r,6)dw\ 

V dr 2 r dr ) 

+ c D (r,0) ( 

d 3 w 

dr 2 dt + 

9) d 2 w \ 
r drdt) 

Me 

= D(r,0) | 

( 1 dw ..d 2 t(F 

) + C£>(r, 6) 

/ 1 d 2 w 
V r drdt 

t n\ w \ 

+ " ( ’’ < VW 


and the piezoceramic patch generated excitation moment is M pe = K^BXpe u {t)* Here c a is 
the viscous (air) damping coefficient, and u(t ) is the voltage applied to the patch. With 
E denoting the Young’s modulus, the spatial variables D = l2 fv ^ v 2 ) ' /y ’ P au< * c ° re P“ 
resent the flexural rigidity, Poisson’s ratio, the mass density per area, and Kelvin- Voigt 
damping for the plate/patch structure. The constant K b is a piezoelectric parameter 
depending on the material piezoelectric properties as well as geometry, and the char- 
acteristic function Xpe is given by Xpe(^) = 1, for r < a pe , Xpe(?') = 0, for for a 

patch of radius a pe . The term V 2 M pe in (20) is an unbounded operator involving Dirac 
delta function derivatives. 

Let fi denote the region occupied by the plate and H pe that of the patch. The 
energy (weak or variational) form of system (20) is 


/ 

Ju 


d 2 w _ , 


+ 



•i£ ii,L ’ + L M ’w dLj+ 

K'b u(t) V 2 r? du> + grjduj 



for a class of test functions g (see [8] for details). 


( 21 ) 


4. PLATE VIBRATION CONTROL 

For the control examples discussed here, we will concentrate on the situation where 
the plate starts with a given initial displacement and velocity and is then allowed to 
vibrate. It is also assumed that there are no external forces applied, i.e. g(t,r,6) = 0. 
The goal in the control problem is to determine a voltage u(t) which, when applied to 
the piezoceramic patch, leads to a significantly reduced level of vibration. For our dis- 
cussions in this note, the control problem formulation ( 1 )-(2) is used and we implement 
the optimal observer of (3)-(8). Robustness of the controller was not considered in our 
first attempt of implementation of the control design. However, it is the subject of our 
current efforts and will be reported on elsewhere. 
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The system describing the dynamics of the plate is infinite dimensional. To approx- 
imate the plate dynamics, a Fourier-Galerkin scheme is used to discretize the infinite 
dimensional system (21). Following the ideas detailed in [8], the plate displacement is 
approximated by 

w N (t,r,9) = ^w^(t)B n (r,d) (22) 

?l— 1 

where { B n (r, #)}(/_, are cubic spline/ Fourier basis functions. The substitution of the 
expansion (22) into (21) yields the 2/V x 2 N matrix system 

y N (t) = A N y N (t) + B N u(t), 

y N m = «!, (2 ' !) 

where y N (t) = (t), • • * , w^t), w ^ (/),***, w^(t)] denotes the 2N x 1 vector contain- 

ing the generalized Fourier coefficients for the approximate displacement and velocity 
(see [8] for details concerning the discretization of the circular plate equation and for- 
mulation of the matrices A ^ and To simplify the notation, the superscript N 

(which is fixed) will be dropped hereafter in this note* The systems in what follows are 
understood to be finite dimensional. 

It has been shown in [5, 6, 12] that the approximation scheme is well defined in 
the sense that solution to the finite dimensional system (23) converges to the solution 
to the original infinite dimensic a 1 system (21). 

For the finite dimensional approximate system, the problem of determining a con- 
trolling voltage can be posed as the problem of finding u(t) which minimizes the cost 
function (2) where y(t) is the solution to (23). From the control design results in 
Section 2, the optimal controller and observer are easily obtained from (3)-(7). 

We report here on a finite dimensional compensator for the approximate system. 
A natural question is whether this compensator will stabilize the infinite dimensional 
system. For bounded input and bounded output operator systems, we refer to [13, 15] 
for detailed discussions on this issue. For the systems with unbounded input and 
output operators, additional results and conditions under which well-posedness and 
convergence are assured can be found in [4, 16, 17$ 18, 20] as well as other references. 

5. NUMERICAL EXAMPLES 

It is known that the PDE-based control design introduced in preceding sections 
requires accurate knowledge of system parameters p, Z), v, c D , c a and K, B < Even 
though material handbooks may provide partial information, the damping coefficients 
are always unknown and the piezoelectric material constant is given only up to certain 
range of values. Before the feedback control law can be designed and implemented, 
significant parameter identification efforts must be carried out. The methodology and 
results for theoretical issues for this parameter identification is reported in [3, 8]. Spe- 
cific parameter identification results using experimental data which were obtained from 
our circular plate were reported in [1]. The same experimental setup was later used in 
control law implementation. The dimension of the aluminum plate and piezoceramic 
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patch are summarized in Table 1. The table also contains “handbook” values for the 
Young’s modulus, Poisson ratio and density of the plate and patch. The estimated 
parameters via fitting model response to the experimental data are summarized in Ta- 
ble 2. As explained in [12], the parameters p, D, c D and v have discontinuities (at the 
patch boundary r — a v ) which must be estimated. 



Plate Properties 

Patch Properties 

Radius 
Thickness 
Young’s modulus 
Density 
Poisson ratio 
Strain coefficient 

a = .2286 (m) 
h = .00127 (m) 

£ = 7.1 x 10 10 (N/m 2 ) 
p = 2700 (kg/m 3 ) 
v = .33 

rad = .01905 (m) 

T - ,0001778(7>i) 

E pe = 6.3 x 10 10 (N/m 2 ) 
p P e = 7600 (kg/m 3 ) 
u pe = .31 

d 3l = 190 x 10- 12 (77(/V ) 


Table 1 . Plate and PZT properties. 



p (kg/m 2 ) 

D(N ■ m) 

cp (N ■■ m ■ s) 
x 10 -4 

u 

C a 

( s • N/m) 

/Cfl 

(N/V) 

beam 

b+P 

beam 

b+P 

beam 

b+P 

beam 

b+P 

Ana. 

3.429 


13.601 




.33 


1 — . 

.013369 

Est. 

3.157 

3.123 

11.017 

11.178 

2.158 

2.210 

.3304 

.3271 

15.566 

.015288 


Table 2. Analytical and estimated values of the physical parameters. 


Using these estimated values of the physical parameters, simulation studies were 
carried out. To closely resemble the experimental setup, we assumed that a single 
point observation, velocity at center, is available. The nonnegative 2N X 2 N matrix Q 
was chosen by taking energy into consideration and weighted as explained in [3], and 
the positive matrix R is just a positive constant which penalizes unrealistically large 
voltages. The matrices Q , R are chosen to be 2 N x 2N and p x p identity matrices, 
respectively, where p is the number of observations ( p = 1 in our simulation study). 

The simulation was carried out in two steps. First, the PDE system with an 
external excitation force and without control (u(t) = 0) was solved for the time period 
of [0,U]- The excitation force was cut off before time U. The solutions at t\ were then 
used as initial conditions (displacement and velocity) in solving the system with control. 
A recorded impact hammer hit was used as the excitation force. The simulation result 
is depicted in Figure 1. In this figure, plot (a) is a time history of uncontrolled versus 
controlled velocities, and plot (b) is the control voltage fed back to the piezoceramic 
patch. The maximum voltage reflected the choice of weights d = 1 and R e = 10 
for the design parameters Q and R respectively. It was observed that it is the ratio 
R e /d which influences the amplitude of the controlling voltage. The sampling time was 
set to 1/12000 Hz, and dimension of the approximation was set to N = 16 under the 
criteria that solution to (23) does not vary significantly if the dimension was larger 
than 16. The time t\ for which the control feedback loop was closed was taken to be 
0.01 seconds. 
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The simulation studies were most encouraging and we subsequently implemented 
the control design in the Acoustic Division at NASA Langley Research Center. A num- 
ber of practical problems arose during the implementation; a significant one involved 
the use of an accelerometer as sensor. To obtain the structural velocity, an approximate 
integrator was used to integrate the acceleration signal. This integrator is defined by 
(see [14]) 

v(t.) + ft v(f.) + ft 2 v(t.) = a(t), (24) 

where v(t) is velocity and a(t) is acceleration. The equation (24) has zero DC gain 
and frequency response close to that of exact integration for u > 6fl. Since the first 
structural mode was approximately 58 Hz, the design parameter D was set to 1 6?r rad/s. 
Equation (24) was implemented by solving the following equivalent first order equation 


v(t) 


’ -n 

ft ' 


«(0 

+ 

1 ' 

. . 


-ft 

0 


v(t) _ 

0 


a(f), 


in which the second variable v was introduced. 

After obtaining the velocity, the estimated state at time t is evaluated by 

y e (t ) = (A - BRT'B ' n - PC'R - 1 C) y c {t ) + PC' RT X v(t) (25) 

where v(t) is the integration of acceleration at the center of the plate. The backward 
Euler method was employed to solve both differential equations (24) and (25). We 
should point out that the last term in (25) is different from the one in (3) which was 
used in simulation. In (3), y(t ) (in y 0 b(t) = Cy(t)) is a vector of generalized Fourier 
coefficients, while in the implementation, the observation y 0 b(t) is an observed variable 
which can be used immediately instead of multiplying by the matrix C. The controlling 
voltage u(t) was obtained through 


u(t) = Ky c {t) (26) 

and then fed back into the system. The matrices (^A — BR~ l B'U — PC'R~ } C^ and K 
in (25) and (26) were calculated offline to reduce computation time in implementation. 

As a preliminary investigation into the feasibility of implementing the feedback 
control scheme, a series of experiments were conducted in which the plate was excited 
with a centered impact hammer strike and the vibrations were recorded both with 
and without controlling voltage being applied. In all cases, the same electronic setup 
was used with the control cases differing only in that the calculated voltage was fed 
back into the system whereas it was simply calculated and stored in the uncontrolled 
case. Representative plots of the plate velocity (integrated from the data recorded by 
the centered accelerometer) for the uncontrolled and controlled cases are given in Fig- 
ure 2. The forces delivered by the hammer impact in the uncontrolled and controlled 
case were nearly identical. By comparing the velocity levels in the two cases, it was 
noted that a significant reduction was obtained in the controlled case with a controlling 
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Velocity (m/sec) 




voltage having a maximum value below 60 volts. Quantitatively, the application of 
the controlling voltage leads to a 47% reduction in maximum velocity levels by time 
t = 0.5 seconds and a reduction of 73% by t — 1 second. We are happy to report 
that these results are typical of those recorded in a series of experiments and hence 
represent the typical reduction in vibration levels that were obtained when the con- 
trolling voltage was fed back into the system in the manner described above. While 
implementation procedures are not yet optimal, these preliminary tests demonstrated 
that vibration levels could be effectively reduced when the PDE-based control scheme 
was implemented. Further experiments are currently underway and will be reported 
on at this conference. 


6. CONCLUDING REMARKS 

In this paper, the experimental implementation of a PDE-based feedback control 
strategy for a flexible structure has been considered. The structure under considera- 
tion was a thin circular plate with a centered piezoceramic patch. When the control 
law was experimentally implemented, the control technique led to a significant reduc- 
tions in the plate vibrations. As we have noted, even more vibration reduction was 
observed in corresponding simulation studies. We are currently conducting research 
to investigate the gap between the simulation and experimental findings. While the 
results of implementation are not yet ideal, they provide a test as to the feasibility of 
implementing a PDE-based method in this manner as well as indicating directions for 
future experimental tests to further improve techniques and methodology. 
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